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Low Energy Neutrino Cross Sections: Comparison of Various Monte 
Carlo Predictions to Experimental Data 

G.P. Zeller'^ 

^Columbia University, Department of Pliysics, New York, NY 10027 

Charged current (CC) and neutral current (NC) low energy neutrino cross section predictions from a variety of 
Monte Carlo generators in present use are compared against existing experimental data. Comparisons are made 
to experimental data on quasi-elastic, resonant and coherent single pion production, multiple pion production, 
single kaon production, and total inclusive cross sections, and are restricted to the case of i^^ scattering off free 
nucleons. 



o 

m 
o 

Oh: 



X 



1. Introduction 

Present atmospheric and accelerator based 
neutrino oscillation experiments operate at low 
neutrino energies {E„ ~ 1 GeV) to access the rel- 
evant regions of oscillation parameter space. As 
such, they require precise knowledge of the cross 
sections for neutrino-nucleon interactions in the 
sub-to-few GeV range. At these energies, neu- 
trino interactions are predominantly quasi-elastic 
(QE) or single pion production processes (Fig- 
ure ^ , which historically have not been as well 
studied as deep inelastic scattering (DIS) reac- 
tions which dominate at higher energies. 




Figure 1. Charged current neutrino cross sec- 
tions as a function of energy (in GeV). Shown 
are the contributions from quasi-elastic (dashed) , 
single pion (dot-dash) and deep inelastic scatter- 
ing (dotted) processes. Figure from Reference pQ. 



production come mainly from bubble chamber, 
spark chamber, and emulsion experiments that 
ran decades ago. Despite relatively poor statis- 
tics and large neutrino flux uncertainties, they 
provide an important and necessary constraint 
on Monte Garlo models in present use. Recent 
neutrino experiments employ a variety of Monte 
Carlo generators to model low energy neutrino 
interaction cross sections. Many share common 
theoretical inputs such as Llewellyn Smith free 
nucleon QE cross sections |7j, Rein and Sehgal- 
based resonance production '16' , along with stan- 
dard DIS formulas and parton distribution func- 
tions. Yet the generators can differ substantially 
in how they implement Fermi gas models, com- 
bine resonance and DIS regions, and treat nu- 
clear and final state effects [T]. This work is an 
attempt to compare several Monte Carlo genera- 
tors to the existing body of low energy neutrino 
cross section data. In particular, three simu- 
lations are considered: v2 NUANCE 2 , NEU- 
GEN 3 , and NUX For the purpose of com- 
parison, the Monte Carlo cross section predic- 
tions in each case have been generated under the 
same set of parameter assumptions _5i, namely: 
ruA = 1.032 GeV, my = 0.84 GeV, gA = ~1.25, 
and sin^ 9\y = 0.233 where applicable. While 
only free nucleon cross section comparisons are 
presented here, detailed comparisons of generated 
Monte Carlo event kinematics are provided in |^. 



Data on quasi-clastic scattering and single pion 
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2. Quasi-Elastic Scattering 

At low energies, CC neutrino hadron interac- 
tions are predominantly quasi-elastic: 

In predicting the cross section for quasi-elastic 
scattering off free nucleons, most Monte Carlos 
commonly employ the Llewellyn Smith formal- 
ism [7] in which the QE differential cross section 
takes the form: 
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where (-1-)— refers to (anti)neutrino scattering, 
Gf is the Fermi coupling constant, is the 
squared four-momentum transfer (Q^ = —q^ > 
0), M is the nucleon mass, E,y is the incident neu- 
trino energy, and (s — m) = AM — Q'^ — . The 
factors A, B, and C are functions of the two vec- 
tor form factors Fi and F2 , the axial vector form 
factor Fa, and the pseudoscalar form factor Fp: 
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where t = Q^/AM^ and m is the muon mass. 
Monte Carlos commonly assume a dipole form for 
the factors Fi, F2, Fa, and Fp: 
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where m^r is the pion mass, fip = 1.793 /xjv and 
fin = —1.913 fiN are the proton and neutron 
anomalous magnetic moments, respectively, and 
the parameters my, tua, and gA are empirical in- 
puts 5]. Departures from this dipole approxima- 
tion which better fit electron scattering data have 
recently been explored and have few-% effects on 
the shape of the predicted cross section . 

Over the years, quasi-elastic processes have 
been studied extensively in bubble chamber ex- 
periments at ANL, BNL, CERN, FNAL, and Ser- 
pukhov. The bulk of this data came from light 
targets and had limited precision due to large 
neutrino flux uncertainties. Figure |2 compares 
the various Monte Carlo predictions to this col- 
lection of QE cross section data. It is no surprise 
that given the same underlying model and input 
parameters, the various Monte Carlo predictions 
agree up to their numerical precision; however, 
note that there is a large spread in the available 
experimental data. 

CC Quasi— Elastic Cross Section 

~g : ■ Serpukov. Belikov, Z. Phys. A320, 625 (1985). Al 
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E„ (GeV) 

Figure 2. QE data compared to various Monte 
Carlo calculations assuming a dipole form for the 
vector and axial vector form factors with my = 
0.84 GeV, niA = 1.032 GeV, and gA = -1.25. 
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3. NC Elastic Scattering 

Neutrinos can also elastically scatter from both 
protons and neutrons in the target material: 

i^^ n — > f ^ n 

Equations I1I4I still apply in describing the neu- 
tral current elastic scattering cross section with 
the exception that in this case the form factors 
include additional coupling factors and a contri- 
bution from strange quarks: 



this ratio from BNL E734 0. While the BNL 
E734 result is quoted over a particular range 
(Table , no restrictions have been placed 
on the Monte Carlo predictions. Also note that 
NEUGEN is 20% larger than the other predic- 
tions. This is simply because NEUGEN assumes 
As = —0.15, which enters the differential cross 
section as (As)^, whereas the other models as- 
sume As = by default \TU\ . 
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where gA = —1.25 and — -|-1(— 1) for proton 
(neutron) scattering. Fl2{Q^) ^re the strange 
vector form factors; the strange axial vector form 
factor is commonly denoted as: 
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where As is the strange quark contribution to the 
nucleon spin. 

Experiments typically publish NC elastic cross 
sections with respect to the CC QE cross section 
to minimize systematics. Table ^ lists a collec- 
tion of experimental measurements of the NC /CC 
ratio, (i/pp v^p)/{Vfj^n tJi~ p)- Figure 01 
shows a comparison between the Monte Carlo cal- 
culations and the most precise measurement of 
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Figure 3. Most precise measurement of the 
NC/CC ratio, [v^p — > v^p)/{v^ n ^ l^^ P)! com- 
pared to predictions from NUANCE, NEUGEN, 
and another neutrino generator, GENEVE I15| . 



Experiment 


Target 


Result 


Q^(GcV^) 


BNL E734 9 


CH2 


0.153 ± 0.018 


0.5 1.0 


BNL CIB 11 


Al 


0.11 ± 0.03 


0.3 ^ 0.9 


Aachen (12| 


Al 


0.10 ± 0.03 


0.2 -> 1.0 


BNL E613 [12] 


CH2 


0.11 ± 0.02 


0.4 -> 0.9 


Gargamelle HH 


CPs Br 


0.12 ± 0.06 


0.3 -> 1.0 



Table 1 

Several measurements of the ratio, {i^p,p 
p) I (i^fj. ^ ~^ p)- Also indicated is the 
interval over which the ratio was measured. 
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4. CC and NC Single Pion Production 

The dominant means of single pion production 
in low energy neutrino interactions arises through 
the excitation of a baryon resonance (N*) which 
then decays to a nucleon-pion final state: 
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N* 
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ttN' 



where A'', N' = n,p. There are seven possible res- 
onant single pion reaction channels, three charged 
current: 
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and four neutral current: 
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Traditionally, Monte Carlos base their theoretical 
calculations of resonant pion production on the 
Rein and Sehgal model [El. While the A(1232) 
is the dominant resonance at these energies, both 
the NUANCE and NEUGEN generators include 
additional higher mass resonant states. In con- 
trast, the NUX model does not yet contain ex- 
plicit resonance production, and is thus why the 
NUX predictions exhibit less agreement with the 
data. Such models are commonly tuned to re- 
produce single pion data, but remain poorly con- 
strained because of the limited availibility and 
large uncertainties in this data. With a few ex- 
ceptions, most of the experimental measurements 
come from deuterium or hydrogen target bubble 
chamber experiments. Figures I4I6I compare the 
Monte Carlo predictions to available CC single 
pion data. The NUANCE and NEUGEN pre- 
dictions include both resonant and nonresonant 
contributions and assume no invariant mass re- 
strictions. Some caution is warranted in draw- 
ing conclusions from these plots as some of the 
higher energy data includes a W < 2.0 GeV in- 
variant mass cut. The data sets with W cuts are 
indicated in each caption. 
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CC Single Pion Production 

■ CERN-WA25, Allasia, Nucl. Phys. B343. 285 (1990), 

• ANL. Barish, Phys, Rev. D19, 2521 (1979). H,, Dj 

• ANL, Rodecky, Phys. Rev. D25. 1161 ( 1 982), Hj, Dj 
o BNL, Kitogoki, Phys. Rev. D34, 2554 (1986). D, 

□ SKAT, Grobosch, Z. Phys. C41, 527 (1989), CF,Br 
A BEBC, Allen, Nucl. Phys. B264, 221 (1986), Hj 

• FNAL, Bell, Phys. Rev. Lett. 41, 1008 (1978), Hj 

• ANL, Campbell, Phys. Rev. Lett. 30, 335 (1973), Hj 
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Figure 4. Measurements of the CC In cross sec- 
tion a{v^p — > /i~p7r+). The Monte Carlos as- 
sume niA = 1.032 GeV and my = 0.84 GeV. 
The data do not include an invariant mass cut 
with the exception of the CERN-WA25, SKAT, 
BEBC, and FNAL measurements which are re- 
ported for W < 2.0 GeV. 
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■ CERN-WA25, Allasia, Nucl. Phys. B343, 285 (1990), D, 
» ANL, Barish, Phys. Rev. 0)4. 2521 (1979), Hj, Dj 
' ANL, Rodecky, Phys. Re/. D25, 1 161 (1982), H„ D, 
o BNL, Kitagaki, Phys. Riv. D34, 2554 (1986), D, 
□ SKAT, Grabosch, Z. Pl^ys. G41, 527 (1989), GFjBr 
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Figure 5. Measurements of the CC In cross sec- 
tion a{v^n — > jj,~nn'^). The data do not in- 
clude an invariant mass cut with the exception 
of the CERN-WA25 and SKAT measurements 
which are reported for W < 2.0 GeV. 
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CC Single Pion Production 
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' '° E.(lfeV) 

Figure 6. CC In cross section cr(i/^ n — > /x~p7r°). 
Same Monte Carlo and data criteria as in Fig- 
ureEl Note: the CERN-WA25 data on this chan- 
nel as reported in the Durham reaction database 
is actually the sum of their measured n — > 



' p tt" and \ 



cross sections [T7] . 



The data on NC single pion cross sections is 
even more limited. Almost all of this data exists 
in the form of NC/CC ratios. Table|21summarizes 
the various measurements, all of which were con- 
ducted using bubble chambers with the exception 
of BNL and CERN PS which both utilized spark 
chambers. In some instances, the experimental 
data can differ by as much as factors of two or 
three. Also listed are the predictions from NU- 
ANCE. In all cases, the NUANCE Monte Carlo 
agrees with at least one measurement. 

Furthermore, data on absolute inclusive NC 
single pion cross sections is extremely sparse. 
While the ANL 12 ft deuterium bubble cham- 
ber experiment ^Sl, reported a cross section for 



the NC Itt channel, n 



v^p-n , the only 



measurements of the remaining NC Itt cross sec- 
tions come from a recent reanalysis of Gargamelle 
propane- freon bubble chamber data [201. Fig- 
ures 171101 compare these measurements to reso- 
nant + nonresonant predictions from NUANCE 
and NEUGEN. Both the data and Monte Carlo 
in all cases have been corrected to free nucleon 
cross sections. 
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NC Single Pion Production 



E : , GGM. Krenz, Nuc. Phys. B135. 45 (1978), CjH.+CFjBr 

O 225 - corrected to o free nucleon cross section 
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Figure 7. NC Itt cross section (j{v^p Vf^pir^). 
Shown are the free nucleon cross section predic- 
tions from NUANCE and NEUGEN with niA = 
1.032 GeV, mv = 0.84 GeV, and sin^ Ow = 0.233. 
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NC Single Pion Production 



£ : , GGM, Krenz. Nuc. Phys. B13b, 45 (1978). CjH.+CFjBr 

U 225 - corrected to a free nucleon cross section 




Figure 8. NC Itt cross section a{v^n v^mr^). 
Same Monte Carlo settings as in Figure [3 
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Table 2 

Measurements of NC/CC single pion cross section ratios. The Gargamelle data has been corrected to 
a free nucleon ratio 22 . Also quoted are the free nucleons cross section predictions from NUANCE 
assuming tua — 1.032 GeV, my = 0.84, and sin^ 9w = 0.2319 in each case. * In their later paper [12], 
Derrick et al. remark that while this result is 1. 6a smaller than their previous result (18) . the neutron 
background in this case was better understood. ** The BNL NC tt" data was later reanalyzed after 
properly taking into account multi-7r backgrounds and found to have a larger fractional cross section |24| . 
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Figure 9. NC In cross section a{i^f_i p ^ v^-i 
Same Monte Carlo settings as in Figured 
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NC Single Pion Production 
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Figure 10. NC cross section (j{v^n v^ipn 
Same Monte Carlo settings as in Figured 
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5. Single Pion Kinematic Comparisons 

So far these comparisons have been restricted 
to the case of free nucleon cross sections; how- 
ever, comparisons can also be made to measured 
kinematic distributions. Such evaluations were 
originally performed in testing NEUGEN gener- 
ator performance and presented at NuIntOl [H]. 
The study is repeated here after including sim- 
ilar NUANCE calculations. Figures I11I12I dis- 
play invariant mass and distributions for the 
three CC single pion channels as measured in the 
BNL 7 foot deuterium bubble chamber . Both 
the NUANCE and NEUGEN Monte Carlo predic- 
tions were generated assuming the BNL flux and 
deuterium target. The Monte Carlo normaliza- 
tion in all plots is determined by the peak of the 
Pfj^p fj," p 7r+ invariant mass distribution in the 
data. Both Monte Carlo models yield comparable 
agreement with the experimental data. 



BNL single pion distributions 



BNL Single Pion Invariant Mass Distributions 
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Figure 11. Invariant mass distributions {W) for 
CC single pion production channels as measured 
in Reference Note the clear A(1232) peak. 
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Figure 12. Four-momentum transfer (Q^) distri- 
butions for CC Itt channels as measured in |26| . 



6. Coherent Pion Production 

In addition to resonance production, neutri- 
nos can also coherently produce single pion final 
states. In this case, the neutrino coherently scat- 
ters from the entire nucleus, transferring negligi- 
ble energy to the target {A). The result is a dis- 
tinctly forward-scattered single pion. Both NC 
and CC coherent pion production processes are 
possible: 



A ~f /i An 



The cross sections for such processes are predicted 
to be small, but have been measured in a variety 
of neutrino experiments. A comprehensive review 
of the experimental data is provided in Refer- 
ence P7) . FiguresHSandlHshow a comparison of 
this data to the NUANCE and NEUGEN predic- 
tions. Both NC and CC data are displayed on the 
same plot after rescaling the NC data, assuming 
o-^^ = 1/2 a^^ I2H1- In addition, data from vari- 
ous targets are corrected to oxygen cross sections 
assuming 

scaling |2H1- 
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CC Coherent Pion Production Cross Section 
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Figure 13. Coherent pion production data. In 
this comparison, NUANCE assumes itia — 1.032 
GeV while NEUGEN uses = 1-0 GeV. 
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Figure 15. Cross section for v^n ^ fjL pn^ tt 
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Figure 14. Same as Figure El except extended 
out to higher energies. 
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Figure 16. Cross section for i/^p ^ /i~p7r+ 7r°. 
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7. Multi-Pion Production 

Because of the complexities involved in final 
state identification and isolation, data on mul- 
tiple pion production is even less abundant and 
less precise than single pion production chan- 
nels. Likewise, such channels exhibit the largest 
spread in the Monte Carlo estimates. Figures IT^ 
1171 compare available data on dipion production 
to NUANCE and NEUGEN. Because deep inelas- 
tic scattering is a large contribution to such cross 
sections, the NUANCE prediction is expected to 
improve with planned DIS model upgrades pg|. 

Multi Pion Production 
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Figure 17. Cross section for i^^p 
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8. Single Kaon Production 

Proton decay modes containing a final state 
kaon, p ^ v , have large branching ratios 
in many SUSY GUT models. Because there is 
a non-zero probability that an atmospheric neu- 
trino interaction can mimic such a proton decay 
signature, estimating these background rates is 
an important component to such searches. The 
following lists some of the contributing strange 
production channels available at low energies: 



v., n 



CC : 



AO 
P 



NC : 

^i^pi^+A" 
>J^pif"AO 



v,,p 



^ n 

fi- K+ E+ 
^l- K+ 



'^p.P- 

v^n - 
v„. n - 



K- E+ 



(10) 



Typically, such reactions have smaller cross sec- 
tions than their single pion counterparts due to 
the kaon mass and because the kaon channels 
are not enhanced by any dominant resonance (in 
contrast to A(1232) decays to single pion final 
states). There are few predictive theoretical mod- 
els for single kaon production PH] and little ex- 
perimental data. Most of the data comes from 
bubble chamber measurements where the strange 
particle decays could be explicitly identified. Fig- 
ure 1181 shows the only two experiments which 
have published cross sections on the dominant as- 
sociated production channel, v^n ^ ^~ K^A^. 
Both bubble chamber measurements were made 
on a deuterium target and based on less than 
30 events combined. To model kaon production, 
NUANCE and NEUGEN employ the same Rein 
and Sehgal-based framework as used to model 
single pion production, including additional reso- 
nance decays other than N* ^ tt N. Note that as 
plotted, neither Monte Carlo includes DIS contri- 
butions to this channel. 
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Figure 18. Measurements of the associated pro- 
duction cross section, cr(z^p n — + /i^ K'^ A^). The 
predictions include resonant contributions only. 
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9. Total CC Cross Section 



C i/„ Total Cross Sections 



Figures I19I20I compare total inclusive CC 
cross section predictions to available experimental 
data. While the total cross section at high energy 
(DIS regime) is known to a few percent, the cross 
section at lower energies is much less precisely 
known. In particular, data measurements in the 
few-GeV range are generally of ^ 10% precision 
and come primarily from experiments which ran 
in the 1970's and early 1980's [21]. At these ener- 
gies, it is especially challenging to model the total 
cross section as there are substantial overlapping 
contributions from QE, resonance, and DIS pro- 
cesses. The Monte Carlo predictions plotted in 
Figure [T^ include all of these contributions, and 
are shown to agree fairly well with each other. 
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Figure 19. Total isoscalar inclusive CC cross sec- 
tion per GeV. The predictions are a sum of all CC 
contributions (i.e., QE, Itt, multi-7r, DIS, etc.). 



10. Conclusions 

This work was an attempt to present a com- 
prehensive comparison between available Monte 
Carlo generators and experimental data. While 
the comparisons are restricted to the case of free 




Figure 20. Same data and Monte Carlo as Fig- 
ure El with inclusion of QE data from Figure El 



nucleon cross sections, the hope is to expand this 
study to include neutrino-nucleus cross sections 
and predictions from additional Monte Carlo gen- 
erators. In the comparisons shown here, reason- 
able agreement between the various Monte Carlo 
simulations and experimental data is observed. 
The agreement is slightly better for CC interac- 
tions than for NC. However, it is certainly true 
that the Monte Carlos agree with the experimen- 
tal measurements to the extent that the data 
agree with themselves. Error bars in the current 
neutrino data set more or less span any observed 
differences in the Monte Carlo predictions. 
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